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Abstract 

We present the first measurements of the nuclear modification factor (Raa) for flavor-separated 
b, c-quark electrons in Au+Au collisions at -Js NN = 200 GeV. The newly installed Silicon Vertex 
Tracker is used to measure the distance of closest approach distributions of electrons at midra- 
pidity (|t7|< 0.35) over the transverse momentum range 1 <p e T < 5 GeV/c. From this, the relative 
fraction of bottom (b) and charm (c) quarks is determined in both the p+p and Au+Au collision 
systems, which form the basis of the measured Raa- In P+P, we observe that a FONLL perturba- 
tive QCD calculation of b — >e/{c — >e+b — ><?) ratio is in good agreement with the data. In Au+Au, 
the data imply a large suppression of fo— >e or a large modification of B meson pj distributions, 
which implies very interesting physics of B mesons in Au+Au collisions. 



1 1. Physics Motivation 

2 In relativistic heavy ion collisions at RHIC, heavy quarks (charm, c, or bottom, b) are ex- 

3 pected to be primarily created from initial hard parton scatterings [ 1 1 and carry information from 

4 the system at an early stage. The interaction between heavy quarks and the medium is sensitive 

5 to the medium dynamics, therefore heavy quarks are suggested as an ideal probe to quantify the 
e properties of the strongly interacting QCD matter. Heavy quark production has been studied by 
7 the PHENIX experiment via the measurement of electrons from semi-leptonic decays of hadrons 
s carrying charm (noted c) or bottom (noted b) quarks. A large suppression and strong elliptic flow 

9 of single electron heavy flavor has been observed in Au+Au collisions at y/s Nf/ = 200 GeV |2). 

10 This suppression is found to be similar to that of light mesons which implies a substantial energy 

11 loss of fast heavy quarks while traversing the medium. The strong flow implies that the same 

12 heavy quarks are in fact sensitive to the pressure gradients driving hydrodynamic flow - giving 

13 new insight into the strongly coupled nature of the QGP fluid at these temperatures. For these 

14 earlier results, PHENIX was not able to distinguish electrons from c and b quarks. In order to 

15 understand medium effects in more detail, it became imperative to directly measure the nuclear 

16 modification, and the flow, of c and b separately. Based on this motivation, in December 2010, 

17 the PHENIX Collaboration opened a new era for measuring heavy flavor at RHIC by installing a 
is new detector called the Silicon Vertex Tracker (VTX). 
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Figure 1: DCA distribution of 2<p e T <2.5 GeV/c electrons measured at -Js^ = 200 GeV in (a) minimum- 
bias Au+Au and (b) p+p collisions. The curves represent the decomposed contribution of charm, bottom, 
photonic background, and other components. The total curve is shown for the fit range used to extract the b 
and c ratio (see text for details). 



The new VTX detector benefits three areas for PHENIX heavy flavor measurements. First, 
by selecting electrons with a distance of closest approach (DCA) to the primary vertex larger 
than ~100/im, the photonic electron background is suppressed by orders of magnitude. This 
suppression results in a clean and robust measurement of heavy flavor production in the single 
electron channel. Secondly, as the lifetime of mesons containing bottom is significantly longer 
than those containing charm, the detailed DCA distribution from the VTX allows to disentangle 
charm from bottom production over a broad pj range. Thirdly, a DCA cut to remove hadrons 
reduces the combinatorial background of Kn such that a direct measurement of D mesons through 
this decay channel will be possible. In this paper, we present the first PHENIX measurements 
using the VTX to measure the DCA of electrons at midrapidity. From this, the yield ratio of 
single electrons from bottom to those from all heavy flavor, as well as the nuclear modification 
factor of charm and bottom separately, are presented using minimum-bias (MB) Au+Au and 
p+p collisions at -Js NN = 200 GeV obtained in the 201 1 and 2012 RHIC runs. 

2. PHENIX Central Detectors 

Electron identification in PHENIX utilizes the two central-arm detectors |3). Until 2010, a 
combination of tracker and electromagnetic calorimeter systems allowed for electron, photon, 
and hadron measurements over the range |?7|<0.35 and A<p-n/2 in azimuth. In 2010, the central 
detector was upgraded to include the VTX ready for Run 201 1. Starting from the beam line, the 
complete detector now comprises the new VTX, which is followed by another tracker system 
made of two sets of drift chambers and a pad chamber (PCI). Outside of the PCI detector is a 
Ring Image Cherenkov detector (RICH), which provides electron/hadron separation from pj of 
a few hundred MeV/c to about 5 GeV/c. The last layer of the spectrometer is an electromagnetic 
calorimeter which is used for photon measurements and electron/hadron separation. 

The VTX was commissioned in Run 201 1 using p+p collision data at y[s — 500 GeV. Com- 
missioning was followed by a period of data collection, which forms the basis of the physics 
measurements presented here, in 2011 (for Au+Au) and 2012 (for p+p) RHIC runs. The VTX 
detector consists of four layers of barrel detectors which cover |?7l<1.2 and almost 2n in azimuth. 
The inner two barrels consist of a silicon pixel device with 50x425 pm pixel size. The outer 
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Figure 2: (a) Fraction of fo— »e to the total (c— >e + >e) as a function of electron p" T for minimum-bias 
Au+Au and p+p at -yJs NN = 200 GeV. Panel (a) shows a comparison between the two systems. These 
extracted fractions were obtained assuming Pythia D and B p T distributions. Au+Au data include no 
uncertainties from modified p T spectra. Panel (b) shows a comparison between p+p and a FONLL |6|. 



47 barrels consist of silicon strip detectors with stereoscopic strips of 80 pm x 3 cm, these devices 

48 achieve an effective pixel size of 80 x lOOOyum. Details of the VTX can be found in Ref. J4). 

49 3. Measurement of the Nuclear Modification Factor for Charm and Bottom 

so The first step in measuring Raa for c and b is to determine the distance of closest approach 

51 (DCA) of the candidate electrons to the primary collision vertex. As electrons from D and B 

52 mesons decays do not originate at the primary collision vertex, they have a large DCA. For ex- 

53 ample, the decay-lengths of D° and B° are 123 and 457 pm respectively. This differentiation 

54 using the DCA is only made possible by the tracking through the VTX. In the present measure- 

55 ments, a precise primary collision vertex — 0-5% central Au+Au collisions at -\Js NN = 200 GeV 

56 has a resolution (cr x , <r y , cr~) of (54±2, 37±2, 68+2) yum — is determined by stand-alone-tracking 

57 using the VTX. These resolutions will be improved by future fine adjustments to the alignment 
se of the detector. 

59 Once candidate electrons are determined from the original PHENIX central arm components, 

eo they are matching to a set of clusters in the VTX. The momentum of the track is measured by 

61 the Drift Chamber. Electron identification is done by RICH and energy momentum matching. 

62 The VTX-portion of the reconstructed candidate is used to determine the DCA to the collision 

63 vertex allows the statistical extraction of flavor identification. Figure [T] shows the raw DCA 

64 distributions of electrons measured at -Js NN = 200 GeV for both minimum bias Au+Au (panel 

65 a) and p+p (b) collisions in the range 2<p e T <2.5 GeV/c. The DCA resolutions are estimated to 
ee be 70 and 138 fan, in Au+Au and p+p respectively. In the present work, the DCA resolution 
67 for p+p was determined relative to the beam center, rather than the collision vertex as used with 
ee Au+Au collisions, which avoids auto-correlations due to the low multiplicity. 

69 The curves presented in Fig.[T]represent the decomposition of the DCA distributions from all 

70 expected signal and background contributions. For each p e T bin, we first define a DCA fit func- 

71 tion, which has contributions from photonic, charm, bottom, K e ^, and a constant-BG, from which 

72 the component yields are obtained. The amount of hadron contamination was determined using 
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Figure 3: Raa of electrons from charm and bottom decays in Au+Au collisions at ^Js NN = 200 GeV. The 
results are compared to published PHENIX n° (7). 



73 the RICH swap method - a standard procedure of PHENIX electron analyses [5|. The DC A 

74 distributions of the photonic components (Dalitz, photon conversion) are Gaussian distributions, 

75 which was confirmed with a full Geant simulation of the PHENIX detector. We additionally 

76 studied the effect of high multiplicities (expected in Au+Au collisions) on the DCA distribution 

77 by embedding simulated electrons and pions in real Au+Au data. This study indicated that the 

78 effect on the DCA due to the high occupancy is small and justifies the use of a pure Gaussian 

79 distribution as the shape of the photonic component. The DCA distribution shape of the c, b, 
so and Kei decay components were determined by simulation. The c and b distribution shapes were 
si generated using Pythia (b—>e and c— »e) decays convoluted with the DCA resolution of the VTX 

82 detector. This was cross checked using a full simulation. The K e i background normalization is 

83 determined from the number of charged hadrons in the same pj bin, the measured on K/h ratio, 

84 the branching ratios of K e i, and a simulation of K e ^ decays. There are additional backgrounds at 

85 large DCA in the data that are caused by conversion electrons in the outer VTX layers. In low 
sb multiplicity events, these outer layer conversions are removed by the requirement of hits in the 
87 first two barrels B0 and B 1 of the VTX. However, in Au+Au collisions, some of these conver- 
se sions are confirmed by random hits in the VTX. Since they are caused by random matching with 
89 unrelated hits in the VTX, their DCA distribution is broad, almost flat. This random background 
go was evaluated using a sideband method which showed that the distribution is consistent with flat 

91 DCA along z in the region 0.75 < |DCA Z | < 1.0 mm. 

92 The total number of inclusive electrons, N'" c , are obtained after subtracting the backgrounds 

93 (hadron contamination, K^, and constant BG) from the total number electrons in the DCA distri- 

94 bution. The number of b and c are obtained using N b = N bK 'xR HQ xF b and JV C = N^' c xR HQ x{\- 

95 F b ), respectively. Here F b is the fraction of b in heavy flavor electrons, F b -b —>e/(c —><?+£> —>e), 

96 which is the only unknown parameter. The Rhq is defined as Rhq = N^ F /N'" c which was es- 

97 timated by two methods: one based on 2004 data of photonic/non-photonic ratio and the other 

98 based on measured conversion tagging probability using 201 1 data. 

99 Figure^ shows the resulting bottom fraction, F b =b — ><?/(c — >e+b — »e), as a function of 

100 electron p e T measured in MB Au+Au and for p+p collisions at -ds m - 200 GeV. The ratio in 

101 p+p collisions is compared to a fixed-order-plus-next-to-leading-log (FONLL) perturbative QCD 
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102 calculation (pink curves) J6) as shown in Fig.JZj}. We observe that FONLL is in good agreement 

103 with current and previous published |5| PHENIX results. Owing to the large uncertainties in the 

104 p+p data the FONLL calculation shape is varied to fit all experimental data points of F h (FONLL 

105 blue solid curve). The blue dashed curves represent lcr uncertainties from the FONLL shape fit 
loo procedure. The F h ratio is used to form Raa for c and b using Eqn.[T| 

RaT-R'aT 6 ^- and R^e =Rb ^elzlAA_ (1) 

ppFit ppFit 

107 where R^T^ e is the nuclear modification factor for single electrons heavy flavor measured by 

108 PHENIX |S). We observe that the nuclear modification of c is less than that for 7r°s (R c A 2 e> R'^ A ), 

109 as shown in Fig. [3] 

no The bottom/charm separation results, shown in Figs|2]and[3] are extracted assuming Pythia 

m D and B pj distributions. These results demonstrate that the Au+Au data are inconsistent with 

112 these input assumptions unless there is also a large suppression of electrons from bottom across 

113 the measured p e T range. This large suppression implies a large change in the parent hadron 

114 pj distributions, which results in changes in the electron DCA distributions. This causes an 

115 additional set of uncertainties in R 1 ^" which are not included in Figure [5] Because the charm 
us fraction dominates at low p T , the R^ is less affected. We are actively working on evaluation 
117 of these uncertainties. These results imply that either a large suppression of b — > e or a large 
us modification of B meson p T distributions, which implies very interesting physics of B mesons in 
no Au+Au collisions. 

120 4. Summary 

121 The DCA distributions of single electrons have been measured by the VTX detector in p+p 

122 and Au+Au collisions at -\Js NN = 200 GeV. By selecting electrons with a DCA to the primary 

123 vertex larger than 100 pm, the photonic electron background was suppressed by orders of mag- 

124 nitude. This suppression results in a clean and robust measurement of heavy flavor production 

125 in the single electron channel. From this, the relative fraction of b and c quarks is determined. 

126 In p+p, the FONLL perturbative QCD calculation of b — > e/(c — » e + b — > e) ratio is in good 

127 agreement with the data. This allows us to measure Raa for single electrons from D and B de- 

128 cays. We observe a higher Raa, i-e- less suppression, for D mesons compared that observed 

129 for 7T°s. While for B, the data imply a large suppression of b — > e or a large modification of B 

130 meson p T distributions. It should be noted that the bottom/charm separation results are extracted 

131 assuming Pythia D and B pj distribution. This leads to additional uncertainties in Raa which are 

132 not included in the present results; PHENIX Collaboration is actively working on evaluation of 

133 the uncertainties. 
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